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Introduction 

The first large gap magnetic suspension systems 
were developed more than 30 years ago as wind 
tunnel magnetic suspension and balance systems 
(MSBS). (See ref. 1 .) More recently, large gap mag- 
netic suspension systems have been developed for vi- 
bration isolation and pointing applications. Differ- 
ences in coil geometries and positioning requirements 
for the vibration isolation application have led to new 
approaches for sensing and controlling the position 
and attitude of the suspended element. In this pa- 
per, the design and performance of an optical posi- 
tion measurement system (OMS) are analyzed. The 
OMS was developed for a six-degree-of- freedom, large 
gap magnetic suspension system (LGMSS) (ref. 2 ) for 
vibration isolation and pointing. 

The LGMSS features several large electromagnets 
in a planar configuration that suspend a cylindrical 
element containing a permanent magnet core (fig. 1 ). 
The element is suspended 91 cm (36 in.) above the 
magnets. To stabilize suspension and control mo- 
tion in six degrees of freedom, information on the 
position and on the attitude of the suspended ele- 
ment is required. An optical sensing system based 
on linear charge- coup led device (CCD) sensors has 
been designed to determine the position and attitude 
of the suspended element in six degrees of freedom 
and to supply this information to the LGMSS con- 
trol system. In the OMS, multiple one-dimensional 
imaging sensors are used to detect small infrared 
light-emitting diode (LED) targets embedded in the 
surface of the suspended element. The position and 
attitude of the element are determined from the mea- 
sured locations of the images in the sensors and from 
transformation equations that relate the coordinates 
of the target images in the reference frames of the 
sensors to the position and orientation of the sus- 
pended element in the laboratory reference frame. 
Rigid body motion has been assumed. 

The position and orientation of the suspended el- 
ement are defined by the alignment of a body fixed 
coordinate system defined by x y and z 5 with re- 
spect to a laboratory fixed coordinate system defined 
by xj , yj, and zj. (See fig. 1 .) Both coordinate sys- 
tems are right-hand Cartesian. The position of the 
element is defined by the x, y, and z location of the 
origin of the body frame in the laboratory frame of 
reference (denoted x cm , y r7n , and z rm ) and the ori- 
entation of the body frame with respect to the lab- 
oratory frame defined by the Euler angles of yaw, 
pitch, and roll and the rotation sequence depicted in 
figure 2 . The magnetic suspension system must po- 
sition the element to xj, yj , and zj — ±0.0254 cm 
(± 0.01 in.) of the initial 0 -in. alignment and to i/>, 6 , 


and (p = ± 0 . 02 °. The LGMSS must meet these spec- 
ifications for a range of element sizes and weights. 
In particular, suspension to this accuracy is to be 
demonstrated for two sizes of cylindrical elements: 
35.56 cm (14 in.) long by 5.08 cm (2.0 in.) in di- 
ameter and 35.56 cm long by 10.16 cm (4.0 in.) in 
diameter. (See fig. 3.) In addition, the LGMSS must 
respond to yaw commands of 1.0° over 360° within a 
10 -sec interval allowed for repositioning and settling. 
No commands can be input for pitch and roll; rather, 
they are constrained to 0 °. 

The LGMSS has a digital controller. Both the ac- 
curacy of the position estimate and the rate at which 
position and attitude information is made available 
to the controller influence the positioning accuracy 
of the magnetic suspension system. An accuracy re- 
quirement of 30 percent of the total error budget for 
the LGMSS was established for the position measure- 
ment system. That is, the position sensing system 
should track the element from xp = 0 ° to 360° with 
an accuracy (i.e., 1 standard deviation) for x cm , yctm 
and z cm of ±0.00762 cm (±0.003 in.) and for xp, 0, 
and <p of ±0.006°. 

The rate at which the sensing system must up- 
date the estimates of position and attitude for the 
controller to stabilize suspension of the suspended 
element is a function of the natural frequencies of 
the magnetic suspension system. The natural fre- 
quencies are functions of the configuration of the coil, 
size and mass distribution of the element and core, 
and orientation of the core magnetization vector. Ini- 
tially, the LGMSS was set at 20 samples per second. 
At the time this sample rate was set, the number 
of coils, the exact coil configuration, the size of the 
suspended element, and the permanent magnet core 
sizes for the LGMSS were unknown. Therefore, the 
final system configuration as designed operates at a 
rate of 40 samples per second with either a small 
or large cylindrical element that contains a small or 
a large permanent magnet core, respectively, which 
has been magnetized in the vertical direction. (See 

fig. 3 .) 

Symbols and Abbreviations 

6 , / coordinate systems of body and 

laboratory fixed reference frame, 
respectively (fig. 1 ) 

F error in position of suspended 

element 

I(n) intensity (fig. 14) 

k time index 



fc, 7 . mu elements of rotation matrices for x 

and y cameras, respectively 

N x , N y window width in pixels of x and 

cameras, respectively 

rfo position vector for point target in 

body reference frame 

r cm position vector of center of body 

reference frame in laboratory 
reference frame 


tj position vector for point target in 

laboratory reference frame 

R position and attitude vector of 

suspended element 

T transformation matrix for rigid 

body (Euler) rotation 

W matrix of weights for target 

positions 

Wj weight for jth target position 

W x , W y weight for centroid estimate in x 

and y cameras, respectively 

x{i, j), y(i, j) centroid location in «th x (or y) 
sensor of target j 

x c , y c , z c camera position in laboratory 
reference frame 


Xfc, yb, Zb coordinates of a point target in 
body reference frame 

x nn , Hem, Zcm coordinates of origin of body 
reference frame in laboratory 
reference frame 


27 , y / , zf 

x 

A 

e 

a 

G 

<t> 

1l> 


coordinates of a point target in 
laboratory reference frame 

centroid of light distribution 

(fig. 14) 

vector of delta changes in position 
and attitude of suspended element 

pitch angle of element, deg 

standard deviation 

calculation of second moment of 
light distribution (fig. 14) 

roll angle of element, deg 

CCD clock signals (fig, 17) 

yaw angle of element, deg 


Abbreviations: 

charge-coupled device 
digital signal processor 
light-emitting diode 

large-gap magnetic suspension 
system 

optical position measurement 
system 

transistor- transistor logic 

Peak signal 

— Standard deviation of background noise 

Design of Optical Measurement System 

Sensors 

The initial design and development of specifica- 
tions for the LGMSS began in 1986. Simultane- 
ously, a feasibility study of the LGMSS was initi- 
ated (ref. 3) and a survey of commercially available 
position measurement systems was conducted. The 
survey results showed that no commercially available 
position measurement system could meet all the re- 
quirements for the LGMSS and that a measurement 
system would have to be developed to meet these 
requirements. Also, optical sensing was determined 
as best for position measurement in the LGMSS be- 
cause of its potential for high-resolution, high-speed, 
and nonintrusive measurements. Laser- based inter- 
ferometric and CCD-based optical sensing systems 
were found to have the greatest reliability, speed, 
and measurement accuracy. Of these two types, the 
CCD-based sensing systems were known to be rel- 
atively impervious to magnetic fields. At the time 
of the survey, a CCD-based position measurement 
system was already in use in a large gap magnetic 
suspension system, the 13-Inch Magnetic Suspen- 
sion and Balance System (MSBS) wind tunnel at the 
Langley Research Center. (See ref. 4.) Preliminary 
findings from the Langley MSBS experiments showed 
that a CCD-based system would provide more reli- 
able and accurate measurements in the LGMSS en- 
vironment than would a laser-based interferometric 
system. 

Two CCD sensor geometries were evaluated for 
use in the LGMSS application: linear CCD arrays 
in which small discrete photosensitive elements are 
arranged in a single line and planar CCD arrays 
in which the photosensitive elements are arranged 
in a grid. The state of the art in processing and 
position estimation was assessed and a close-range 
photogrammetric approach was selected for track- 
ing the location of the suspended element. In this 
approach, multiple cameras are used in conjunction 


CCD 

DSP 

LED 

LGMSS 

OMS 

TTL 

SNR 
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with multiple-point targets and triangulation to de- 
termine the locations of the targets in an external 
frame of reference. This selection was based upon the 
fact that more sophisticated algorithms which relied 
upon full images of the suspended element or pat- 
tern matching were too computationally intensive, 
given the capability of processors at that time, to 
meet the rate requirement of 20 samples per second. 
Therefore, the CCD sensor geometries were evalu- 
ated based on the number of cameras required to 
meet the accuracy requirements for point tracking, 
and the speed or time required to acquire data from 
the array and to track the different targets. 

To estimate the number and spatial resolution of 
cameras required to meet the LGMSS accuracy speci- 
fications, a photograrnmetric sensing system was sim- 
ulated using the simultaneous triangulation and re- 
section system (STARS) computer program. (See 
ref. 5.) Cameras were located above; the suspended 
element and angled as close to 45° as possible with 
respect to the vertical. (The errors in the estimates 
°f xj, yj , and zj are a minimum for a 45° pointing 
angle. See ref. 6.) For a fixed camera resolution and 
a fixed uncertainty in the camera locations in the 
external frame of reference, the accuracy of detect- 
ing a point target within the field of view improves 
approximately as the square root of the number of 
cameras. Assuming that the location of a target 
image could be resolved to 0.1 pixel in a camera, 
6 two-dimensional CCD array cameras with 1024 by 
1024 pixels would be required to yield a 1 a un- 
certainty in the position of a point target in the lab- 
oratory reference frame of 0.001016 cm (0.0004 in.) 
in xj and yj and 0.001524 cm (0.0006 in.) in zj. 

At the time of the evaluation, the state of the 
art for planar CCD arrays was a device with 512 
by 512 pixels. Thus, more than 20 two-dimensional 
CCD cameras would be needed to achieve the same 
accuracy in point tracking as was predicted for a 
camera with 1024 by 1024 pixels. However, state-of- 
the-art linear arrays featured 4096 pixels. Based on 
the assumption that two linear CCD arrays oriented 
orthogonally would have a resolution equivalent to a 
two-dimensional array, only four linear CCD cameras 
with 4096 pixels would be needed. (This assumption, 
in hindsight, was a poor one. System calibration 
of linear CCD arrays is more difficult and the un- 
certainties in the camera parameters can be larger 
than were modeled in the STARS simulation.) 

At least three targets are required to solve the 
transformation equations. Then, by assuming equal 
uncertainty in the locations of targets and equal 
moment arms, the analyst can improve the estimates 
of position and attitude roughly as the square root of 


the number of targets. Thus, for a la uncertainty in 
a target position of 0.001016 cm in xj and yj and of 
0.001524 cm in zj, at least six targets are needed to 
meet the yaw and pitch accuracy requirements. (To 
meet the accuracy requirement for roll for the smaller 
element, more than 10 targets are required, assuming 
the same la uncertainty in target position.) 

In the case of the two-dimensional array, all the 
targets could be imaged simultaneously and, there- 
fore, the only consideration was whether the image 
data could be digitized and stored fast enough to 
allow time for processing. However, the linear ar- 
ray would require the targets to be multiplexed in 
time for target identification. Therefore, the total 
time would be the product of the number of targets 
and the time to acquire a single target. In CCD 
arrays, the signals generated by each photosensitive 
element (or pixel) are read out in a serial fashion. 
As a result, for a fixed data rate, the greater the 
number of pixels in the array, the greater the read- 
out time. The maximum data rate for CCD arrays 
at the time of the evaluation was about 20 MHz. 
This data rate corresponded to the conversion rates 
for the fastest analog-to-digital converters (ADC’s) 
(above 8- bit resolution). At these rates, an array 
with 512 by 512 pixels required at least 13 msec to 
be read out. The time required to read out a linear 
array of 4096 pixels at 20 MHz was 205 //sec. There- 
fore, assuming an integration time of about 2 msec, 
the time to image and store the data from six targets 
was about 15 msec for the two-dimensional array and 
13 msec for the linear array. 

The choice of CCD array geometry was ultimately 
based on the greater one-dimensional spatial resolu- 
tion of linear CCD arrays and the fact that the de- 
sign of the electronics for data storage and processing 
was much simpler for the linear array. Thus, a linear 
CCD array with 2048 pixels was selected as the basic 
sensor element for the optical measurement system 
(OMS). This CCD array had the maximum pixel den- 
sity for linear CCD arrays at that time. The system 
was designed for 16 cameras with 8 targets embedded 
in the surface of the suspended element. This design 
ensured that at least 6 targets could each be resolved 
to 0.1 pixel in at least 12 cameras for the full range 
of motion of the suspended element. 

Targets 

The decision was made early in the OMS design 
phase to use light-emitting diodes (LED’s) as active 
targets rather than passive reflective targets. Passive 
targets do not require electronics or power on the 
suspended element; however, they must be large to 
return a strong signal. Passive targets require a 
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powerful source of illumination; also, when a linear 
CCD detector is used, one passive target is not 
readily distinguishable from another. Further, the 
large target centroid location is more difficult to 
determine from oblique camera angles. 

Developing active targets for the OMS proved 
challenging for several reasons. If the targets were 
to be detected from different camera perspectives, 
then the light from the targets would have to radiate 
in a nearly hemispherical pattern. This approach 
would require high light output from the LED’s to 
obtain a high camera signal-to-noise ratio (SNR) for 
a short integration period (less than 3 msec). To 
obtain that high a light output, the LED’s would 
have to be driven at the maximum allowable current 
levels. High current levels tended to overheat the 
LED’s; also, because the LED’s were to be powered 
with batteries, the run time decreased. 

Figure 4 shows a schematic diagram and photo- 
graphs of the LED targets developed for the OMS. 
A single LED chip that measures 0.05 cm (0.02 in.) 
across the diagonal is soldered into a reflector with an 
outside diameter measuring 0.11938 cm (0.047 in.). 
The base of the reflector shown in figure 4 is a heat 
sink. Initial tests indicated that the heat sink was 
required to prevent destruction of the LED when it 
was operating at the high current pulses used in this 
application because early tests had disclosed that 
nearly 50 percent of the light was attenuated by a 
diffuser. During tracking, the LED’s are pulsed on for 
2.083 msec at 200 mA. The LED’s are top emitting 
with a wavelength of 880 nm. Furthermore, no lens 
or defuser is used on the LED. Figure 5 show's the 
target geometry and a plot of the normalized target 
light output as a function of angle. 

The LED’s are multiplexed in time for target 
identification. Power and driver electronics for the 
LED’s are located inside the cylinder. (See fig. 6.) 
Two rechargeable batteries are used to power the 
targets and electronics. The batteries are rated for 
250 mA- hours of operation. Tests on a single battery 
showed that battery lifetime will depend upon the 
method of charge and discharge. Curves of output 
voltage over time are shown in figure 7 for different 
discharge rates. 

Processors 

The design of the processor evolved over a 5-year 
period. This evolution was driven by the require- 
ments for high speed and high accuracy. Maximum 
accuracy and speed required that the centroid esti- 
mate be calculated using floating point arithmetic. 
The total time to transfer the data from the cam- 


eras to a single central processor, to execute the cen- 
troid algorithm 16 times in sequence, and to solve the 
transformation equations exceeded 50 msec. As a re- 
sult, a decision was made to design an ADC board 
and a digital signal processor (DSP) board for each 
camera. The digitized data from a camera would be 
stored and the centroid of the one-dimensional target 
image would be calculated on the DSP board. The 
16 centroids would be transferred to a central proces- 
sor where the x, y , and 2 positions of the targets and 
the position and attitude of the suspended element 
would be calculated. 

In 1987, standard DSP’s and microprocessors did 
not have built-in floating point capability. Thus, 
additional floating point processors were originally 
included in the design of the camera processors. 
Advancements in technology, which occurred around 
1988, led to the use of high-speed DSP’s with built-in 
floating point capability. 

The design of the central processor evolved as 
well. Originally, the calculations of the target posi- 
tions and the position and attitude of the suspended 
element were to be performed in a high-speed array 
processor that was external to the OMS. However, 
communications between the OMS and external pro- 
cessor proved too slow. The basic processor was fast 
enough to perform the computations in the required 
time, but the data transfer times between the sen- 
sor DSP’s and the external processor were too great 
to meet the 20-sample-per-second requirement. As 
a result, DSP 17 (MATH DSP) was added to the 
OMS backplane and used to calculate the target po- 
sitions and the position and attitude of the suspended 
element. 

The software for the processor is structured so 
that data acquisition and computations are done 
both in parallel as well as pipelined. A diagram il- 
lustrating the structure of the software is shown in 
figure 8. Each frame is divided into 12 intervals, the 
intervals being equal to the camera integration peri- 
ods. During interval 1, the first target is pulsed on 
and imaged in parallel in the 16 cameras. During 
interval 2, the data are digitized, stored in random 
access memory (RAM) on each camera DSP board, 
and the calculations of the centroids are begun in 
the camera DSP’s as the second target is pulsed on 
and imaged. At the end of interval 2, the centroids 
are transferred to the MATH DSP. During interval 3, 
the third target is imaged, the centroids of the second 
target are calculated in the camera DSP’s, and the 
x , y, and 2 position of the first target is calculated 
in the MATH DSP. Intervals 4 8 are identical to in- 
terval 3. During intervals 9-12, the targets are not 
illuminated. The calculations of the centroids and 
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the x , ?/, and z position estimates of the eighth tar- 
get are calculated during intervals 9 and 10, respec- 
tively. During the last two intervals, the position and 
attitude of the suspended element are calculated. 

Optical Measurement System 

The OMS consists of eight sensing units stationed 
symmetrically about the center of the suspension 
magnets and 8 ft above the top of the magnets. (See 
fig. 9.) Each sensing unit has two linear charge- 
coupled device (CCD) cameras, denoted x and y , 
oriented in orthogonal directions. (See fig. 10.) The 
sensing units are attached to a support structure that 
is constructed of wood and fiberglass beams with 
aluminum mounts. The mounts allow each sensing 
unit to be oriented in two independent angular di- 
rections. The distance between each sensing unit 
and the origin of the laboratory reference frame is 
approximately 183 cm (72 in.). 

Camera components include a Thompson-CSF 
linear CCD array sensor with 2048 photosensitive el- 
ements, a 5.08- by 5.08-cm single-element cylindri- 
cal lens, an infrared filter, driver electronics for the 
CCD array, and thermoelectric cooler and control cir- 
cuitry. The cooler and control circuitry is required to 
keep electronic components on the CCD driver board 
from overheating at the clock speeds used in this 
application. 

The cylindrical lens images the light emitted by 
an LED target as a one- dimensional (line) image on 
the CCD array. Charges are generated in each photo- 
diode element in response to light striking the array 
and they accumulate during the sensor integration 
time. Following the integration time, the accumu- 
lated charge is read out as an analog video signal. 
That is, the voltage signal is output for each photo- 
diode element which is proportional to the integrated 
light intensity falling on that element (pixel). The 
sensor field of view is approximately 60 by 60 cm. 
Thus, an LED target covers roughly 0.002 of the to- 
tal field of view. The width of the lens point spread 
function (PSF) is about the width of an individual 
pixel and both contribute to blurring of the target 
image. The target is small relative to the total field 
of view so the image of a target diode when located 
in the center of the field of view appears as a nearly 
Gaussian distribution of light along the array. (See 
fig. 11.) 

The CCD array is not equally sensitive to all 
wavelengths of light; rather, it has a peak sensitivity 
in the infrared around 775 nm. To improve the SNR, 
an infrared filter is mounted in front of the lens to 
limit the background light reaching the sensor. The 
cutoff for the filter is 830 nm. 


An infrared laser is located on the camera support 
structure directly above the suspended element. The 
laser is used to time the flashing of the LED targets. 
A small phototransistor, approximately 0.127 cm 
(0.050 in.) in diameter, is embedded in the center 
of the top surface of the cylindrical element. The 
target LED s are turned on and off in response to 
the illumination of this phototransistor. The flashing 
of the timing laser is controlled through the optical 
sensing system controller. A timing diagram is shown 
in figure 12. 

A 20 megasample per second, 10-bit, ADC digi- 
tizes the signal from a camera and stores the digital 
values in RAM located on the camera DSP board. 
The centroid of the one-dimensional digital image 
is computed using a Texas Instruments TMS320C30 
DSP. In addition, an SNR is determined for each 
measurement. The SNR provides information to de- 
termine weighting factors, which are carried along 
with the position information. In the MATH DSP, 
the position information is multiplied by the weight- 
ing factors to determine an estimate of the position 
and attitude. As the SNR for any given measurement 
approaches 0, the less the measurement contributes 
to the position and attitude estimate. Each x and y 
centroid estimate and weighting factor is stored as a 
32-bit floating point number. 

Figure 13 shows a block diagram of signal flow 
and processing for the OMS. The function of each 
component shown in the diagram of figure 13 is as 
follows. A DSP is used for each camera to control 
the camera background adjustment, integration time, 
and signal readout. The DSP’s are also used to 
perform calculations such as computing the centroid 
of each image formed by an LED target, the SNR 
of the signal, and the background level, as well as 
to perform an estimate of the total error for the 
corresponding centroid estimate. The MATH DSP 
is used to calculate the matrix elements, execute 
the time interpolation algorithm, and perform the 
tracking algorithm to determine the six position and 
attitude estimates. A 16-bit microprocessor is used 
to control the flow of instructions from the DSP’s 
to and from a personal computer and to control 
communication with the LGMSS controller. The 
personal computer is used to send instructions to the 
optical sensing system and to display system status 
information. 

Calculating Position and Attitude 

The x and y locations of the projected images of 
each of the 8 target LED’s in each of the 16 cam- 
eras are determined by calculating the location of 
the centroid of the one-dimensional image (or light 
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distribution) along the respective x and y CCD ar- 
rays. The six parameters, Xcm, Vcm > z cm, an< ^ 

0, which give the position and attitude of the sus- 
pended element at one point in time, are determined 
by solving two sets of linear equations: the first re- 
lates the measured x and y locations of the target 
LED images in the 16 cameras to the xy, t/y, and zj 
coordinates of each target LED; the second relates a 
change in the xy, yj, and zj coordinates of all eight 
target LED’s to a change in the six (position and at- 
titude) parameters. To correct for element motion 
between the image acquisitions of each target LED, 
a Lagrangian time interpolation algorithm is used. 
The interpolation is performed over the xj, y/, and 
zj coordinates of each target using the current es- 
timates and the estimates from the three previous 
frames. The coordinates of each target are predicted 
at the beginning of the current frame. The inter- 
polated xy, t/y, and zj of the target LED s are used 
to predict the position and attitude of the suspended 
element. 

Centroiding Algorithm 

A flowchart of the centroiding algorithm is shown 
in figure 14 in which three searches are conducted 
for the peak pixel location. First, a coarse search of 
every third pixel is made for the peak pixel location, 
(A search of every pixel would take too long.) A 
coarse search for the peak is performed each time 
because targets will move in and out of the sensor 
field of view as the position and attitude of the 
element changes. Searching the entire array for 
the peak requires less software than tracking the 
centroid estimates and determining when a target 
has departed or entered the field of view. If the peak 
location is determined to be within 50 pixels of either 
edge of the array, then an error flag is set and the 
centroid is not computed. 

Second, if the peak lies within the central region 
of the array, then a search of every pixel is made 
for the peak location ±10 pixels on either side of 
the initial peak location. This new peak location is 
stored and used as a pointer into the array. 

Each pixel has a slightly different gain and zero 
point. Before the centroid is computed, the values 
of light intensity are corrected for pixel nonuniform 
responsivity and the background light level is sub- 
tracted for the 81 pixels centered about the peak. 

Third, a search for the peak location is performed 
over the corrected light intensity values. This “true” 
peak location then replaces the old pointer. 


When the peak location is established, the width 
of the target image distribution is determined by 
calculating the second moment of the distribution. 
The width of the distribution is determined because 
the amplitude of the target image will change as 
the target moves within the camera field of view. 
(See fig. 15.) As discussed in reference 7, when 
the signal amplitude and shape vary, a centroid 
calculated over a window that is a fixed size will cause 
a significant error either due to truncation of the 
signal when the image distribution is large relative to 
the window width or to noise when the image width 
is narrow relative to the window width. To increase 
the dynamic range of the imaging system, the second 
moment of the distribution is calculated before the 
centroid and an optimal window; that is, a window 
is determined that is matched in width to the target 
image distribution. The x and y centroid locations 
are subsequently determined by 




n=Peak + 

E n/ hj( n ) 

n=Peak — ^ L 
n=Peak + ^F 

E r hj(n) 

7i=Peak — AgF 


(la) 




V 

m= Peak + 

E mlij(m) 

ra=Peak - ^ 
m=Peak 

E 

rn — Peak — ^ 


(lb) 


where i denotes the particular sensing unit; j denotes 
the target; n and m are the pixel numbers in the x 
and y cameras, respectively; /(n) and I(jn) are the 
digital values of the light intensity for pixels n and 
rn; peak is the pixel location of peak signal; and N x 
and N y are the half-widths of the window in the x 
and y camera, respectively. 


Tracking Algorithm 

The position and attitude of the suspended ele- 
ment are defined by the parameters x cm , y C m , z cm » 
jp, 6 , and <j>. The position vector r b j of diode j in 
the body frame with coordinates x/ ; j, y^j, and Zfy is 
related to the position vector r jj of diode j in the 
fixed frame with coordinates Xfj. yjj , and zjj by 


Y fj — r cm ± 


( 2 ) 


6 


where T, the rotation matrix, is 


T = 


cos(0) cos(0) 
cos(0) sin(v^) 
—sin ( 0 ) 


— cos(0) sin(0) + sin(0) sin(0) cos('ip) 
cos(0) cos(0) + sin(0) sin(0) sin(^) 
sin(0) cos(0) 


sin(0) sin('0) + cos(0) sin(0) cos(^) 
— sin(0) cos(^) + cos(0) sin(0) sin(0) 
cos(0) cos(0) 


Similarly, the x t j and y 2 y coordinates of diode 
j in the frame of sensing unit i are related to the 
coordinates of diode j in the laboratory reference 
frame by the colinearity equations (ref. 8) 


x p — fx 

mi 1 {xjj - x ( - x ) + m v2 {yfj - y ( ix ) + m i; 


"‘21 (x/j ~ + rn 22 (yfj - !/‘ x ) + W 25 

P/j-V J 



( 3 a 


■ yp - fy 


k\ { {x/j -x< iy ) + k x 2 (y fj - y ( : y ) + ku(z fj - z^) 


k2\ (Xfj + *22 (Vfj - l fly) 4 - htiZfj - t y / J 

(3b) 


where y?., and z\ x arc the x, y , and 2 positions 
of the perspective center of camera x (or camera y for 
subscript y) in the fixed frame; } x and f y are the focal 
lengths of the lenses in cameras x and y, respectively; 
x pi and y pi are the principal points of the ith x and y 
cameras; and and kjj are elements of the rotation 
matrices for each camera and are functions of the 
camera pointing angles. The equations are functions 
of the locations of the ith camera in the fixed frame, 
the focal lengths of the cameras, and the fixed-frame 
coordinates of the jth LED target. 

From equations ( 3 a) and ( 3 b), the xy, yy, and 
zy locations of each target are determined using 
the method of weighted least squares. Lagrangian 
interpolation is performed on the current estimates of 
x fi Vf, a nd zj and the three previous time estimates 
to yield a prediction of the locations of each target 
at the beginning of the sample interval. The inter- 
polated x j , yy, and zy locations of the eight targets 
are used to estimate position and attitude. 

The position and attitude of the suspended el- 
ement are determined by solving for the transfor- 
mation from the body reference frame to the labo- 
ratory reference frame, where rigid body motion is 
assumed. The equations are solved using a Gauss- 
Newton method. The solution for sample time k cor- 
responding to the beginning of the fcth frame of data 
is the solution to 


min 1 , , 0 

^j;^\W k (6F k A k + F k )\q (4) 


where 

F k = (17 - r cm - Tr b ) k 

and 6 F k denotes the Jacobian of F k and W k is the 
weighting matrix composed of eight weighting terms 
onc for each target position estimate. The so- 
lution to equation ( 4 ) is obtained using Householder 
QR factorization. (See ref. 9 .) The position and at- 
titude for sample time k are then determined by 

R-fc = R-fc-i + &k 

where R is a vector with components x rrn , y, w , z nn , 
i£\ 8 , (p , and the position and attitude of the model 
at k or k — 1; and A is the solution to equation ( 4 ) 
at k. 


Weighting Factors 

Weighting factors are determined for each x and y 
centroid estimate as well as for each target position 
estimate xj, yj, and zj. Weighting the measure- 
ments can improve the estimate of target position 
when the weights arc inversely related to the error 
in the measurement. Also, weighting the measure- 
ment provides a simple means of eliminating bad data 
(e.g., target not in the field of view, signal too low, 
etc.) When the calculated weight falls below a preset 
threshold, the weight is set to 0 and the measurement 
does not contribute to the estimates of position and 
attitude. 

The calculated weighting factors for the centroid 
estimates are the inverse of the total predicted error 
multiplied by a constant. That is, 

WxiiJ) = ^§J) (5a) 




c 

£y(ij) 


( 5 b) 


where W x (i,j) and W y (i,j) denote the weighting 
factors; £ x {i, j) and e y {i,j) denote the predicted error 
for the ith x or y sensor of the jth target; and C 
denotes the constant. The predicted error £ x (i,j ) or 
e y (i,j) is calculated as the sum of three terms: the 
first term is dependent upon the SNR, the second 
term is the residual error after the cameras were 
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calibrated to determine the optical distortion, and 
the third term is a function of the proximity of the 
current target image to the target images acquired 
during the previous integration periods of the current 
frame. The third error term, the proximity term, is 
calculated to account for any error from the residual 
charge on the CCD sensor. The dependence of the 
error on SNR, residual charge, and residual errors 
after camera calibrations is described in the section, 
“Performance Analysis.” 

When weighting factors are determined for each 
centroid estimate, they can be determined for each 
target position estimate at k. These weighting fac- 
tors are the same for xj , yj , and Zf for a given target 
and sample time and are calculated from the centroid 
weighting factors as 

1 8 

Wj = ±^lW x {iJ) + Wy{iJ)] (6) 


Calibration Techniques 

Calibrations are performed to correct for errors 
and determine the (unknown) locations of each cam- 
era in the laboratory reference frame. Two princi- 
pal calibrations are necessary — calibration of each 
camera and calibration of the overall sensing system. 
Each camera is calibrated to determine both the re- 
sponse of the CCD sensor to the location of light 
along the array and the correction for optical distor- 
tion. Calibration of the overall system is performed 
to determine the exterior orientation parameters of 
the cameras; that is, the locations and pointing an- 
gles of the cameras with respect to the laboratory 
reference frame. 

Camera Calibrations 

Variations in the predicted location of the cen- 
troid as a function of position along the CCD ar- 
ray will result due to nonuniformities in pixel sensi- 
tivity in the detector and optical distortions in the 
projected target image. Calibrations are performed 
to correct for these errors. To determine the pixel 
responsivity for a CCD array, a (nearly) spatially 
uniform direct current (dc) light source is positioned 
over the end of a camera with the lens in place. The 
average output at each pixel is obtained for nine in- 
tegration times ranging from 1 to 9 msec. The aver- 
age output at each integration time is computed for 
100 scans at that setting. A line is fitted to the out- 
put of each pixel as a function of integration time. 
The slope of the line and the Y intercept are calcu- 
lated for each pixel as the averages from 30 trials. 


The slope corresponds to the responsivity and the 
Y intercept corresponds to the voltage output for no 
light on the sensor. The response of each pixel is 
normalized to the peak response along the CCD ar- 
ray. The inverse of the normalized number for each 
pixel and the Y intercept are stored in look-up tables 
and used to correct subsequent scans. Representative 
graphs of sensor nonuniform responsivity correction 
and offset are in figure 16, which shows two distinct 
lines. The distinct lines occur because, as illustrated 
in figure 17, the CCD array has two separate read- 
out lines. One line is connected to and reads out the 
even-numbered pixels, the other line reads out the 
odd-numbered pixels. Because the relative gain and 
offset of the analog signals from these two lines are 
not the same, the average normalized response and 
offset for the two sides are different. 

Figure 18 shows the calibration setup to deter- 
mine the correction for optical distortion. Each cam- 
era is mounted in turn in a test fixture 167.6 cm 
(66 in.) above a single LED. The LED is mounted in 
the center of a plate, which is fixed in position (with 
dowels) to the top of a computer-controlled linear 
X-Y stage. The camera test fixture consists of a 
flat plate that is leveled to approximately 0.004 cm 
per meter. A plumb bob, positioned in the cam- 
era test fixture, is used to align the camera with the 
LED. To perform a camera calibration, the plumb 
bob is removed and a camera is positioned on top 
of the plate so that the CCD array is aligned with 
the X stage. The LED target is stepped across from 
-25.4 to 25.4 cm in X and the centroid of the pro- 
jected LED image is calculated at 301 equally spaced 
positions along the X-axis. The X- and T-axis stage 
displacements are measured simultaneously with a 
laser interferometer. 

Figure 19(a) shows a graph of the centroid esti- 
mate as a function of the stage position and a third- 
order curve fit to the data for sensor 3. Figure 19(b) 
shows the optical distortion calculated from the data 
in figure 19(a). From these measurements, a third- 
order polynomial curve fit is determined for each 
camera and relates the true position of the projected 
target image along the array (in inches) to the calcu- 
lated centroid position. The coefficients of the poly- 
nomial are stored for each sensor and used to correct 
the centroid estimate for optical distortion in sub- 
sequent camera measurements. 

System Calibration 

A system calibration is performed to determine 
positions and angles, the principal point, and the 
focal length of each of the 16 cameras. The cal- 
ibration LED and the computer-controlled linear 
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stages are shown in figures 20(a) and 20(b). A 
single diode is mounted to a plate attached to the 
computer-controlled X-Y stage and aligned with 
dowels. The LED is scanned in a grid measuring 
10.16 by 10.16 cm. A spacer 2.54 crn (1 in.) thick 
and with matching dowel holes is then mounted be- 
tween the stage platform and the LED plate and the 
LED is again scanned in a 10.16- by 10.16-cm grid. 
Then, a second spacer of the same thickness as the 
first is added and the procedure is repeated. In this 
way a grid of points is generated for a volume mea- 
suring 10.16 by 10.16 by 5.08 cm. 

The origin of the laboratory reference frame is 
defined as coincident with the location of the LED 
when the stage is in the 0 reference position (0 in. 
and 0°) and a single spacer has been added. The 
x and y locations of the control points used during 
system calibration are taken to be the laser x and 
y measurements. The 2 locations of the control 
points are taken to be the thickness of each spacer. 
The corresponding camera centroid estimates of the 
control points are the mean centroid estimates of 
100 scans at each x and y stage position. The 
camera parameters are the solution to the generalized 
least-squares problem that minimize the error in the 
laser measurements as well as the camera centroid 
estimates subject to constraint equations that fix the 
distance and angles between an X and Y sensor pair. 

Performance Analysis 

Five dominant factors affect the estimates of sus- 
pended element position and attitude: 

1. SNR’s of the cameras 

2. Residual errors after camera calibrations 

3. Proximity of one target to another within the 
same frame of data 

4. Error in the exterior orientation parameters of the 
cameras either due to residual error after system 
calibration or to motion of the camera support 
structure 

5. Motion of the model 

Two simulations have been constructed to study the 
sensitivity of the position and attitude estimates 
to these five factors. One simulation models the 
digital signal of the target images and is used to 
determine the effect of the SNR and of sampling on 
the accuracy of the centroid estimate. The other 
simulation models the sensing geometry and tracking 
algorithms used in the OMS and is used to study the 
sensitivities of the position and attitude estimates to 
the five factors. The effect of each of the five factors 


on the final estimate of model position is discussed 
next. 

Camera Signal-to-Noise Ratio 

Experiments on the OMS have determined the 
relationship between the SNR and the error in the 
centroid estimate. Figure 21 shows a plot of the total 
error, which is the sum of the standard deviation 
of the centroid estimate and the bias error, as a 
function of the SNR. The figure shows that the 
measured error is approximately inversely related 
to the SNR. Thus, as the SNR increases, the error 
decreases and asymptotically approaches a value of 
about 0.04 pixel. 

A simulation has been developed to investigate 
the influence of the SNR and the corresponding 
change in the weighting factor of each target on 
the accuracy of the position and attitude estimates. 
In the simulation, the SNR for a given target was 
determined by calculating the difference in distance 
between a target and a camera, by normalizing with 
a reference distance, then by correcting for the falloff 
in light as a function of the difference between the 
camera pointing angle and the target pointing angle. 
The positions and angles of each target relative to 
the body reference frame are shown for the smaller 
element in figure 22. No other sources of error were 
assumed in the simulation. Figures 23(a) and 23(b) 
show the predicted errors in x ani yem, z nil , 0 , 0 , 
and 0 as a function of 0 for ifi = 0 ° to 90°; nominal 
values of x cr ni Vcm-i Zcm ~~ 6 in . 5 0 and 0 — 0 °: and 
the maximum SNR = 200. Similar graphs are shown 
in figures 24(a) and 24(b) but with 9 and 0 = 5°, 
and x crn i ycim and Zcrn — 0 in. In figures 25(a) 
and 25(b) are graphs of the error in the position and 
attitude of the cylinder, respectively, as a function of 
the standard deviation of the centroid estimate. All 
16 cameras were assumed to have the same sensitivity 
and noise cKar act eristics. 

Residual Errors After Camera Calibration 

The residual errors after calibration for optical 
distortion were found to be less than 0.01 pixel 
for all cameras. Similarly, provided a sensor has 
not shifted out of calibration, the residual error 
due to camera nonuniform responsivity has been 
found to be less than 0.01 pixel. Table I lists the 
standard deviations of the 16 camera gain and offset 
corrections for responsivity. The magnitude of the 
standard deviations is a function of the difference in 
the gain and offset of the two sensor readout lines. 

The OMS simulation showed that the error in the 
position and attitude estimate due to residual camera 
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errors of 0.01 pixel or less is too small relative to the 
error due to the SNR to be distinguished from the 
SNR error. This situation will be the case provided 
the camera voltages do not drift and the calibration 
for responsivity does not change. However, if a sensor 
shifts out of calibration, that is, if the analog signals 
from the two lines shift relative to one another and 
the sensor is not recalibrated, then a residual system- 
atic error larger than 0.01 pixel can be introduced in 
the centroid estimate. This systematic error is spa- 
tially periodic over the array with a period of two 
pixels. (See fig. 26.) The magnitude of this error will 
be a function of how far out of calibration the sensor 
has drifted. In addition to degrading camera perfor- 
mance, a systematic camera error of this magnitude 
can introduce errors in system calibration. Both er- 
rors will in turn introduce errors in the position and 
attitude estimate. Therefore, the cameras and as- 
sociated electronics must be housed in a thermally 
stable environment to minimize voltage drift. 

Proximity Error 

As described earlier, charges are generated in each 
photodiode element of a CCD array in response to 
light striking that element and accumulating dur- 
ing the sensor integration period. The way in which 
charge is accumulated and transferred depends upon 
the design of the CCD sensor and can have a pro- 
found effect upon the output signal. Figure 27 shows 
a schematic diagram of the charge transfer process 
for the CCD sensor used in the OMS. After the inte- 
gration period, charge that has accumulated at the 
photo sites (photodiode elements) is transferred in 
parallel to serial shift registers at the start of a read 
cycle. After this transfer, the next integration pe- 
riod begins and shifting of the charge to the output 
is initiated. 

The efficiency with which charge is transferred 
along (and out) of the shift register is characterized 
by the charge transfer efficiency (CTE) of the CCD 
array. Typical values of CTE for modern CCD array 
designs are greater than 99.999 percent. Thus, very 
little signal loss occurs in the shift register. More im- 
portant in this case is the efficiency of charge trans- 
fer from the photo site to the shift register. Tests of 
the OMS revealed a low photo site transfer efficiency 
(PTE) of the CCD’s used in this application. As a 
result, not all of the charge is transferred from the 
photo site to the shift register during a single read 
cycle. Rather, some residual charge remains at the 
photo site that, along with the charge accumulated 
during subsequent integration cycles, is transferred 
during subsequent read cycles. Figure 28 shows a 
plot of the peak output signals for an LED target 


image for eight consecutive read cycles in which the 
LED was on only for the integration period preced- 
ing the first read cycle. The plot of figure 28 is like 
that of a capacitive discharge curve. Similar plots 
have been made for different pixels along the same 
CCD array. The curve of figure 28 was found to vary 
little from pixel to pixel; rather, it showed the de- 
cay of residual charge as a function of the number 
of read cycles that characterize the CCD array as a 
whole. For the CCD arrays used in the OMS, several 
read cycles are required to clear the array of residual 
charge. 

The manifestation of photo site residual charge is 
an afterimage or “ghost image” of a target. To clear 
the array of residual charge, multiple read cycles are 
executed at the end of the frame. The charge is then 
allowed to build to the background level with no tar- 
get LED’s illuminated. However, the time between 
the sampling of each target within a frame is in- 
sufficient to clear the CCD array. As a result, after- 
images of the eight target LED’s accumulate during 
any single frame. A digital filter has been incorpo- 
rated in each sensor DSP to filter out these after- 
images. (See fig. 29.) But, as figure 29 shows, the 
filtering operation is not perfect because, when two 
target images are in proximity, an error will be in- 
troduced in the centroid estimate. The magnitude of 
this proximity error will depend upon the amplitude 
of the afterimages and on the distances between the 
afterimages and the current target image. 

An empirically derived function (fig. 30) is used 
to estimate the magnitude of the proximity error 
between the current target and any other target 
acquired within the frame. A total proximity error 
is determined to be the sum of these individual 
estimates and is added to estimates of the SNR 
and residual errors. Then, a weighting factor is 
determined for the total error estimate. 

The OMS simulation was run to analyze the 
effect of the proximity error on the accuracies of 
the position arid attitude estimates. Figures 31(a) 
and 31(b) show the error in position and attitude of 
the element as a function of r/;. The nominal position 
of the element is 0 in. with 0 angular displacement 
and the sensors are assumed to have a maximum 
uncertainty due to noise of 0.03 pixel. No other errors 
have been assumed. 

Errors in Exterior Orientation Parameters 

Errors in the exterior orientation parameters of 
the cameras are dominant sources of error in posi- 
tion and attitude estimates. Such errors will intro- 
duce both bias and random errors in the position 
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and attitude estimate of the suspended element for 
the following reasons. Sixteen camera parameters per 
sensing unit are solved during resection and triangu- 
lation. As a result, within the uncertainty bounds 
of the measurements (i.e., the x, y, and z control 
point locations and camera centroid estimates), sev- 
eral possible solutions can exist. Constraint equa- 
tions have been introduced to limit the solutions to 
those that are compatible with the known camera 
geometry. Nevertheless, the camera parameters can 
vary, and the resulting cost function remains small. 

To study the sensitivity of the position and at- 
titude estimate to errors in the camera parameters, 
mean zero Gaussian random variables were added to 
the camera parameters and the simulation was run 
more than 3000 times to determine the mean and 
standard deviation of the resulting position and atti- 
tude estimates. Table II lists the standard deviations 
of the position and attitude estimates for different 
magnitudes of error added to the x, y, and 2 posi- 
tions and the pointing angles of the cameras. The 
element was in the 0° position and for all runs and 
was assumed to be perfectly stationary. Table III 
lists the standard deviations of the position and at- 
titude estimates for different magnitudes of error 
added to the x, y, and 2 positions of the cameras and 
an angular error in the pointing angles of 0.0005°; 
for table III, the element was in the 0° position and 
motion was included. A trajectory was assumed for 
levitation generated by the LGMSS simulation. 

Tables II and III show that the accuracies of the 
position and attitude estimates are sensitive to er- 
rors in the camera parameters. For the la value 
of position and attitude to remain within the spec- 
ification for the OMS, the standard deviations of 
the errors in the camera parameters must be small 
(i.e., <0.0003 cm (0.0001 in.) in x r , y r , and z r and 
<0.001° for the pointing angles). 

Camera motion during target acquisition will in- 
troduce additional errors in the position and attitude 
estimate of the suspended element. Unlike the resid- 
ual error after system calibration, which introduces 
both a bias and a random error in the position es- 
timate, the error introduced by camera motion will 
vary in time and will be a function of the modes 
of vibration of the camera support structure. Fi- 
nite element analyses have been performed to pre- 
dict motion of the support structure. Accelerometer 
measurements were made of the floor of the building 
that will house the LGMSS and the power spectral 
density of the disturbance input was calculated. This 
information, together with the finite element model, 
yielded a maximum predicted displacement along the 


2 direction of approximately ±0.0004 cm (0.0002 in.) 
for the first mode. 

Model Motion 

Because the target LED s are turned on sequen- 
tially, the position of each target LED is determined 
at a different point in time. In the tracking algo- 
rithm, the suspended element is assumed stationary 
as the positions of the eight target LED’s are ac- 
quired. To correct for target motion during signal 
acquisition, a time interpolation algorithm is used. 
The location of each LED within a frame is inter- 
polated in time based on the four previous frames; 
then, the positions of all eight LED’s are estimated 
at the time the first LED within the frame was turned 
on. Without a correction for motion of the suspended 
element between acquisition of each target, the error 
in the position and attitude estimate grows rapidly 
as a function of integration time and rate of element 
displacement. Even with a correction, the error in 
the position and attitude estimate grows as a func- 
tion of integration time and frame time. 

As an example of the effect of model motion on 
the error in the position and attitude estimate, ta- 
ble IV shows the standard deviations of predicted 
errors in the trajectories of the element based on a 
model of the plant and control system for the mag- 
netic suspension system under three conditions: sus- 
pension about the 0° position, release by the me- 
chanical injection device, and a commanded yaw 
displacement at 1.0° per second. Figures 32(a) (f) 
show the input trajectory for the second case, which 
involved release by the mechanical injection device. 
For the results depicted in table IV, the time be- 
tween frames is 25 msec and the integration time is 
2 msec. A minimum la uncertainty of 0.03 pixel 
was assumed for the centroid estimate; and un- 
certainties of 0.0003 cm and 0.0001° (la) were as- 
sumed in the x, y. and 2 positions of the cameras. 
A comparison of the results of tables II and III, 
which were obtained assuming no motion for the 
model, with the results of table IV, in which inter- 
polation was used, shows that motion of the model 
increases the errors in the estimates of position and 
attitude by about a factor of 2. 

Experimental Results 

Several experiments have been run on the OMS to 
determine system accuracy and the effects of different 
calibration parameters and errors on the final system 
performance. Three principal experiments were run 
to assess overall system accuracy. In the first experi- 
ment, the average xj , yj, and zj positions of a single 
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LED are measured at 100 locations in a 10 by 10 grid 
measuring 10.16 by 10.16 cm. 

The position measurements of the first experi- 
ment are shown in figures 33 and 34. Figures 33(a) 
and 33(b) show the differences between the xj or yj 
OMS-predicted LED coordinates and the laser sys- 
tem plots of the Xf or yj values. The values of 
z predicted by the OMS over the whole grid are 
shown plotted in figure 33(c). The standard devi- 
ations of the Xf, y/ , and zj estimates ranged from 
<0.0005 cm (<0.0002 in.) to 0.001 cm (0.0004 in.) 
over the grid of points. The maximum standard devi- 
ations of Xf, yj, and zj predicted by the OMS sim- 
ulation over the same grid for a minimum la un- 
certainty of 0.03 pixel for the centroid estimate and 
a la uncertainty of 0.0003 cm in the camera param- 
eters are 0.001 cm (0.0005 in.) for Xf and yj and 
0.002 cm (0.0008 in.) for Zf. 

In the second experiment, a half-model of the 
cylinder is mounted on top of a rotary stage that is 
affixed to the x and y linear stages. (See fig. 34.) The 
yaw, pitch, and roll axes are aligned with a laser; and 
the x, y, and 2 positions of the model are determined 
from the x and y laser interferometer measurements 
and the distance between the center of rotation and 
the origin of the coordinate system. The body 
coordinates of the targets are taken to be the xj, yj, 
and zj OMS-predicted target locations of the model 
aligned at \p, 0 , and <j> = 0° and at a displacement 
of x and y = 0 in. for the linear stages. The body 
coordinate system is assumed to be coincident with 
the laboratory coordinate system and, after the Xf , 
yj, and Zf locations of the targets are determined, 
these values are read into the tracking algorithm as 
the y\>, and 2& target locations. The model is 
then moved =L1.27 cm (0.5 in.) in x and y and 
measurements are taken over a 10 by 10 grid to 
determine OMS accuracy for small displacements. 

The results of the second experiment are shown in 
figures 35 and 36. The error in the estimates of Xcm, 
y cin . and z C m versus the difference between the laser 
measurements and the OMS estimates is plotted in 
figure 35. Figure 36 shows the standard deviations of 
the estimates. The standard deviations of the OMS- 
predicted position and attitude estimates using a 
minimum la uncertainty of 0.03 pixel in the centroid 
estimate and a la uncertainty of 0.0003 cm in the 
camera parameters are x cm = 0.0013 cm (0.0005 in.), 
Vcm = 0.0015 cm (0.0006 in.), 2 cm = 0.0025 cm 
(0.0001 in.), ip = 0.002°, 9 = 0.004°, and <p = 0.048°. 

In the third experiment, the repeatability of the 
OMS is determined. The element is aligned at x , y , 
and 2 = 0 in. and ip, 9 , and <p = 0°. The rotary stage 


is set to 0° and the OMS output is noted in all six 
degrees of freedom. The element is rotated in 10° 
increments to ip = 180° (on the rotary stage) and 
the position and attitude are recorded for each yaw 
location. The element is rotated back to ip = 0° 
and the position and attitude are recorded. This 
procedure is repeated 20 times and the positions and 
attitudes are recorded for each yaw location. The 
repeatability of the OMS was found to be better than 
x, y, and 2 = 0.0025 cm (0.001 in.); ^ and 0 » 0.005°; 
and <p ~ 0.01°. 

Concluding Remarks 

An optical measurement system (OMS) has been 
designed and tested for the large gap magnetic sus- 
pension system (LGMSS) that is being used at Lang- 
ley Research Center. The LGMSS is to be used 
to test control laws for magnetic suspension sys- 
tems for aerospace applications of vibration isola- 
tion and pointing. The hardware design, signal flow, 
and details of the tracking algorithm and calibra- 
tion techniques have been examined; also, system re- 
quirements were reviewed and a performance analysis 
done based on simulations and experiments. 

The OMS is based on point target-tracking tech- 
niques and uses 16 linear charge-coupled device 
(CCD) cameras to detect the locations of small in- 
frared light-emitting diodes (LED’s) on a cylindri- 
cal element that is equipped with a permanent mag- 
net core. The location of the projected image of an 
LED in a camera is determined as the position of 
the centroid of the light distribution falling on the 
CCD array. High-speed analog-to-digital converters 
and digital signal processors are used to digitize the 
video signal and compute the centroid location from 
the digitized signal. Based on the locations of the 
projected images of the LED’s in the cameras, the 
sensing system estimates the position and attitude 
of the suspended element and supplies this informa- 
tion to the LGMSS controller. 

The OMS has been calibrated and experiments 
have been conducted to evaluate system accuracy. 
Experiments designed to determine the accuracy of 
point tracking revealed that the system can deter- 
mine the x, y , and 2 location of a target within 
0.001 cm (0.0005 in.) within a volume of 10.16 by 
10.16 by 5.08 cm (4 by 4 by 2 in.). Based on a simula- 
tion model of the OMS, camera signal-to-noise ratios 
consistent with measured values on the actual sys- 
tem, and an uncertainty in the camera parameters 
of 0.0003 cm (0.0001 in.), the predicted standard de- 
viations of the axial estimates of points within the 
volume were 0.001 cm (0.0005 in.). Results of initial 
experiments on a model of the suspended element 
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showed the system in the 0 reference position (0 in. 
and 0°) to be accurate within 0.0025 cm (0.001 in.) 
for the origin of the body in the laboratory reference 
frame, within 0.005° for pitch and yaw, and within 

0. 1. for roll. These same experiments revealed that 
the position and attitude estimate is sensitive to the 
position and attitude of the suspended element and 
that the accuracy will degrade for small pitch or roll 
angles and clement translations >1 cm. The simula- 
tion and the experiment agreed closely in predicting 
the position and attitude of the suspended element. 

Additional tests should be conducted on the sys- 
tem to understand the effect of calibration and test 
techniques on the measured accuracies of position 
and attitude. In addition, a full analysis should be 
made to determine the effect of motion of the sus- 
pended element. 

NASA Langley Research Center 
Hampton, VA 23681-0001 
March 11. 1994 
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Table I. Standard Deviations of Gain and Offset Corrections for All 16 Cameras 


Gain 

0.78 

.43 

.27 

.32 

.27 

.42 

.27 

.99 

2.5 

1.7 

.51 

.36 

1.9 

.28 

.33 

.36 

Offset 

4.2 

2.3 

2.3 

2.6 

1.3 

5.1 

1.5 

2.2 

12.0 

5.4 

2.8 

2.2 

2.8 

1.6 

1.4 

1.9 


Table II. Standard Deviations Predicted by Simulation for Different 
Values of Assumed Uncertainty in Camera Positions 

[Camera resolution was assumed to be 0.03 pixel for peak SNR and cylinder was stationary ] 


No motion for 


Standard deviation of error for 



A positions, cm 

A angles, deg 

Kerri) Cm 

2/cm; cni 

^crfi j cm 

'll), deg 

8, deg 

4>, deg 

0.0003 

0.0001 

0.0005 

0.0008 

0.0003 

0.001 

0.002 

0.02 

.0003 

.0005 

.0008 

.0010 

.0003 

.002 

.003 

.03 

.0005 

.0001 

.0008 

.0008 

.0003 

.001 

.002 

.02 

.0005 

.0005 

.0008 

.0010 

.0003 

.002 

.003 

.02 

.0003 

.001 

.0013 

.0015 

.0008 

.003 

.004 

.04 


Table III. Standard Deviations Predicted by Simulation for Different Values of Assumed Uncertainty 
in x, y, and 2 Positions of Cameras Given That Cylinder is Undergoing Motion 


Simulated input trajectory was steady-state suspension; angular uncertainty was 0.0005° 
for all three cases; camera resolution was assumed to be 0.03 pixel for peak SNR 


Standard deviation of uncertainties 
in position of cameras, cm 

Standard deviation of error for — 

%cm 1 cm 

yem 1 cm 

Zcm 1 cm 

0, deg 

0, deg 

<i>, deg 

0.0003 

0.0013 

0.0013 

0.0005 

0.004 

0.004 

0.05 

.0004 

.0013 

.0013 

.0018 

.008 

.015 

.07 

.0005 

.0015 

.0015 

.0020 

.009 

.018 

.09 


Table IV. Standard Deviations for Position and Attitude Variables Predicted 
by Simulation for Different Simulated Trajectories of Suspended Element 

Camera resolution was assumed to be 0.03 pixel for peak SNR; standard deviation of error in x, y , and z 
positions of cameras was assumed to be 0.0003 cm; angular uncertainty was 0.0005° in all cases 


Conditions 

Id uncertainty for — 

cm 

y, cm 

z, cm 

deg 

8 , deg 

<P » deg 

Levitation 

0.0013 

0.0013 

0.0005 

0.004 

0.004 

0.05 

Launch 

.0014 

.0014 

.0004 

.004 

.004 

.05 

Commanded yaw displacement of 1° 

.0014 

.0015 

.0004 

.002 

.004 

.05 
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Zf, z 
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Figure 3. Orientations of core magnetization vectors and sizes of cylindrical elements. 
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a Variable time required to process image. 

Figure 8. Software structure and timing used during tracking. 



Figure 9. LGMSS and OMS. 
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Figure 11. Digital output of camera 3 expanded to show signal distribution of LED target. Image was acquired 
with camera mounted 64 in. directly above LED target. 
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Figure 12. Target illumination sequence and time allocated for computing position and attitude of suspended 
element. 
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Figure 14. Centroiding algorithm. 
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Digital signal amplitude, counts 




Normalized gain correction 



(a) Offset correction for sensor 5. 


2000 



Pixels 


(b) Normalized gain correction for sensor 5. 
Figure 16. Representative sensor correction and offset 
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Figure 17. Camera analog circuitry and two output signals. 
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Figure 18. Experimental setup with overhead camera, LED. and computer-controlled linear stage. 
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(a) Calibration setup. 

Figure 20. Experimental apparatus. 
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140 
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usor 7. 




o 10 20 30 40 50 60 70 80 90 


Yaw angle, deg 

(a) Predicted error as function of yaw angle (maximum centroid accuracy assumed to be 0.03 pixel). 



0 10 20 30 40 50 60 70 80 90 


Yaw angle, deg 

(b) Predicted error as function of yaw angle. 

Figure 23. Nominal element position of x cm, ycm, and z^n = 0 in.; 6 and <£ = 0°. 


33 



angular error, deg 






(a) Predicted error as function of standard deviation of error in centroid estimate. 



Yaw angle, deg 

(b) Predicted error as function of standard deviation of error in centroid estimate. 

Figure 25. Experiment with all cameras assumed to have same sensitivity and noise. Nominal position of 
element is in all degrees of freedom. 
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Figure 27. CCD array charge transfer circuitry and timing. 
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Figure 30. Empirically derived function used to estimate proximity error of current target sample as function 
of difference in distance between current target sample and target after image and difference in acquisition 
time between current sample and target after image. 


40 




0 10 20 30 40 50 60 70 80 90 


Yaw angle, deg 

(a) Predicted error in position estimate of cylinder as function of yaw angle. 



0 10 20 30 40 50 60 70 80 90 


Yaw angle, deg 

(b) Predicted error in attitude of cylinder as function of yaw angle. 

Figure 31. Assumed minimum centroid error of 0.03 pixel; 0.01 pixel assumed for residual optics error; proximity 
error included. 
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Time, sec 


(a) Predicted x crn trajectory of cylinder during launch by mechanical injection device. Element assumed 
displaced in x initially, then moves to origin of laboratory reference frame under influence of control system. 



Time, sec 

(b) Predicted y cm trajectory of element during launch. 
Figure 32. Input trajectory after release by injection device. 
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Time, sec 


(c) Predicted z rm trajectory of element during launch. 




Time, sec 

(d) Predicted yaw trajectory of element during launch. 
Figure 32. Continued. 
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(a) Difference between average x an measurement from OMS and average laser measurement over ±1.25-cm 
displacement in x and y for 10 by 10 grid. 



(b) Difference between average y cm measurement from OMS and average laser measurement over ±1.25-cm 
displacement in x and y for 10 by 10 grid. 

Figure 35. Three-dimensional plot of half- model experiment. 
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(c) Average z cm measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 grid. 



(d) Average yaw measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 grid. 

Figure 35. Continued. 
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(e) Average pitch measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 grid. 



(f) Average roll measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 grid. 

Figure 35. Concluded. 
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(b) Standard deviation of y (m measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 
grid. 

Figure 36. Three-dimensional plots showing standard deviations for half-model experiment. 
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(d) Standard deviation of yaw measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 
grid. 

Figure 36. Continued. 
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(e) Standard deviation of pitch measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 
grid. 



(f) Standard deviation of roll measurement from OMS over ±1.25-cm displacement in x and y for 10 by 10 
grid. 

Figure 36. Concluded. 
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